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The physics of pulses in gamma-ray bursts: emission 
processes, temporal profiles and time lags. 
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ABSTRACT 

We present a simple, semi-analytical model to explain GRB temporal and spectral 
properties in the context of the internal shock model. Each individual pulse in the 
temporal profiles is produced by the deceleration of fast moving material by a com- 
paratively slower layer within a relativistic wind. The spectral evolution of synthetic 
pulses is first obtained with standard equipartition assumptions to estimate the post- 
shock magnetic field and electron Lorentz factor. We get E p oc t~ s with S = 7/2 
which is much steeper than the observed slopes 5 bs ^ 1-5. We therefore consider the 
possibility that the equipartition parameters depend on the shock strength and post- 
shock density. We then get a much better agreement with the observations and our 
synthetic pulses satisfy both the hardness-intensity and hardness-ffuence correlations. 
We also compute time lags between profiles in different energy channels and we find 
that they decrease with increa sing hardness. We finally compare our predicted time 
lag - luminosity relation to the lNorris et alJ 1)2000(1 result obtained from 6 bursts with 
known redshift. 
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1 INTRODUCTION 

Cosmic gamma-ray burst (hereafter GRBs) exhibit a 
great diversity of duration and profiles. The distribution 
of durations is clearly bimodal with two peaks at about 
0.2 and 20 seconds. GRB light curves are highly variable 
but can often be interpreted in terms of a succession of 
eleme ntary pulses which possibly overlap dNorris et alJ 
1996). These pulses appear as the building blocks of the 
profiles and understanding their physical origin would 
certainly represent a clue for a better description of 
the whole GRB phenomenon. The pulse temporal evo- 
lution has often been described by a fast rise followed 
by an exponen t ial de cay (the so-called FRED shape; see 
iFishman et alJ dl994l) ) but other mathematic al behaviors 
such as stretched expone ntials, gaussian i Norris et alJ 
Il996l) or power-law decays llRvde fc Svenssori l200oir iiave 
been also proposed. Spectral hardness decreases during 
pulse decay and two relations between the temporal and 
spectral propertie s , the H IC (hardness-intensity correlation; 
iGolenetskii et a i] Jl983 |)) and the HFC (hardness-fluence 
correlation: iLiang^^Kargatisl lll996T) ) appear to be satisfied 
by a substantial fraction of GRB pulses during the decay 
phase. Pulse profiles peak earlier in higher energy bands 
and the corresponding time lags between different energy 
channels correlate to pulse hardness and peak luminosity 
jNorris et alJ l200(j) . These observational results must be 



reproduced by the models and may help to discriminate 
among different possibilities. 

Two distinct mechanisms have been proposed to explain 
the origin of pulses in GRBs. In the external shock model 
they are formed when a relativistic shell ejected by the 
central engine is deceler ated by the circumstellar material 
jMeszaros fc Reeslll99l) . An homogeneous medium leads 
to a single pulse but an irregular, clumpy environment can 
produce a c omplex profile if a large n umber of small clouds 
are presgritj|DCTmer_fc_^R^ In the internal shock 

model jRees fc MeszarosHl994D" the central engine generates 
a relativistic flow with a highly non uniform distribution 
of the Lorentz factor and the pulses are made by collisions 
between rapid and slower parts of the flow. In the two 
scenarios the variability of the profiles has a very different 
interpretation. In one case it provides a "tomography" of 
the burst environment while in the second it reveals the 
activity of the central engine. 

In this paper we consider in some details the mechanism 
of pulse formation by internal shocks. Three characteristic 
time scales may be relevant during pulse evolution: the time 
irad required to radiate the energy dissipated in shocks; the 
dynamical time ta yn , i-e. the time taken by internal shocks 
to travel throughout the flow and the angular spreading 
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time tang corresponding to the delay in arrival time of 
photons emitted from a spherical shell. A short radiative 
time t r ad <C tdyn,tang appears to be mandatory to avoid 
adiabatic losses and maintain a sufficient efficiency. This 
condition is satisfied by the synchrotron process which is 
the most commonly invoked radiation mechanism in GRBs. 
If the thickness of colliding shells is small compared to 
their initial separation, td yn <C tang and the pulse temporal 
evolution is fixed by geometry; conversely if the source 
produces a continuous wind rather than a series of discrete, 
well separated shells, td yn > t ang and hydrodynamical effects 
control the pulse shape. 

Pulse evolution has been stud ied extensively when i t 
is dominated by geome try (see e.g. iFenimore et al.l |l996); 
iKobavashi et all 1199711 but discre pancies between model 
predi ctions and the observations ISoderberg fc Fenimord 
l200lT) have cast some doubt about the validity of the in- 
ternal shock model. Our purpose is to see if the situation 
can be improved when the hydrodynamical point of view is 
adopted. We first summarize in Sect. 2 some basic informa- 
tions regarding pulse temporal and spectral evolution. We 
then develop in Sect. 3 a simple model where pulses are 
formed when a fast moving wind is decelerated by a com- 
paratively slower shell. Spectral evolution is considered in 
Sect. 4 where constraints are obtained on the GRB radiation 
mechanism. Temporal profiles computed with our model are 
presented in Sect. 5 and time lags are discussed in Sect. 6. 
Sect. 7 is the conclusion. 



2 TEMPORAL AND SPECTRAL EVOLUTION 
DURING PULSE DECAY 

We consider a pulse characterized by a photon flux N(t) 
in the energy range (Ei, E2), a peak energy E p (t) of the 
E 2 JV(E, t) spectrum and a photon fluence defined by 

r-E, 



d$ N (t) 
dt 



= N(t) = 



Af(E,t)dE 



Ei 



The HIC and the HFC are then given by 
E p (t) oc N(t) s 
and 

E p (t) oc e - a *N(*) ) 



(1) 



(2) 



(3) 



where a is an exponential decay c onstant. For pulse s satis- 
fying both the HIC and the HFC. lRvde fc SvenssorJ (1200(f) 
have shown that the photon flux and the peak energy follow 
simple power laws during the decay phase 
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l+t/7 



N(t) = 
and 

Ep{t) = (T+W ' 

where t — corresponds to the maximum of N(t). 



ip,0 



(4) 



(5) 



iRvde fc Svensso'nl (|2002f l performed a detailed analysis 
of the decay behavior of a sample of 25 long and bright 
pulses to check whether it was indeed described by eqs 2] 
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Figure 1. Derivative of ip(E p ) for three Band functions with 
/3 = -2.5 and a = -2/3 (dotted line) a = -1 (full line) and 
a = —1.5 (dashed line). The two vertical lines limit the BATSE 
spectral range. The average slope during pulse decay typically lies 
between and 1. 

and They found that to account for the temporal and 
spectral evolution of all the pulses, eq. |3]had to be replaced 
by the more general expression 

iVo 



N(t) = 



(1 + t/r)" ' 
If n 7^ 1, the S i ndices appearing in eqs[ 



(6) 



]and|n]are dif- 
ferent and following Rvd e fc Svenssonl J2002T) we then write 
E p (t) = £ p ,o/(l + t/r) 5 * with <5* = nS. IRvde fc Svenssonl 
( 2002) found that the distribution of n in their sample was 
sharply peaked at n — 1 with however a secondary bump at 
n < 3. The values of J* were all smaller than 1.5 for the n = 1 
pulses but could reach 3.5 when n ~ 3 . The distribution of 
5 was narrower with 0.5 < S < 1 in most of the sample. 

Once the decay behavior of N(t) and E p (t) has been 
specified, it becomes possible to obtain the evolution of the 
bolometric energy flux i*E(f) since 

/•oo poo 

Fv(t)= / M(E, t)EAE — E p / N(x,t)xdx , (7) 
J Jo 

where x = E/E p . We suppose that the temporal and spec- 
tral behavior can be separated in N(x, i) : 

M{x,t) = A(t)B(x) , (8) 

B(x) representing the spectrum shape. The photon flux in 
the energy range (Ei, E2) is then given by 

J Ex E p (t)ipa 

so that 

iVo£p,o ^0 



(9) 



F E {t) = N(t)E p (t) 



9ti 



<p(E p ) (l + t/r)i+«* V {E V ) 
with 

p(E p )= / B{x)dx 
Je 1 /e p 

and 



(10) 



(11) 
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Figure 2. Solution for 7(t) corresponding to eg. 1211 with 70 
0.25. The dashed line is the approximation given by ea. 1231 



Figure 3. Dissipated power in units of 2 1C for our pulse model 
with 70 = 0.25. The dashed line corresponds to ea. 1241 while the 
full line takes into account angular spreading. 



ipo = / xB{x)dx . (12) 
'0 

The derivative of <p{E p ) has been represented in Fig. 
for the BATSE spect ral range (20 1000 keV), using a stan- 
dard Band function feand et al.l ll993) with low and high 
energy indices a = —2/3, —1 or —1.5 and /3 = —2.5. At 
low (resp. high) E p , tp(E p ) is given by a simple power- 
law E p '' 3+1 ) (resp. Ep fa+1 ^) but for intermediate values 
(Ei < Ey < E2) which are rep resentative of the decay phase 
in the lRvde fc Svenssonl i2002l) sample, <fi(E p ) does not have 
a simple analytical form. Assuming that it can still be ap- 
proximated by a power-law, <fi(E p ) oc E p where £ is 
a weighted average of the low and high energy spectral in- 
dices (— 2 <C,$ — 1) the bolometric energy flux also follows 
a power-law 



F E (t) oc 



with 



{1 + t/rY 



(13) 



(14) 



(15) 



e = n + (2 + C) <5* • 

The slope of the HIC is then given by 

s= ** = 1 

n e/8* - (2 + C) ' 

If the temporal and spectral evolution during pulse 
decay i s due to geometri cal effects alone e = 3 and 
8* = 1 jGranot et al.lll999l) which leads to 8 = With 
—2 < £ < — 1 the resulting value 0.3 < 8 < 0.5 lies below what 
is fou nd in most observed pulses ^Soderberg fc Fenimord 
|2C 



Geometrical effects govern pulse evolution if the shell 
thickness is small compared to their initial separation. But 
if a continuous outflow emerges from the central engine the 
hydrodynamical time scale can play a dominant role du- 
ring pulse decay. We have then developed a simple model 
to check whether a better agreement can be found with the 
observations when the hydrodynamical aspect of the flow is 
taken into account. 



3 A SIMPLE PULSE MODEL 

We consider a relativistic wind where a slow shell of mass Mo 
and Lorentz factor To decelerates a more rapid part of the 
flow characterized by a constant mass flux M (in the source 
frame) and Lorentz factor Ti > To- We do not solve the true 
hydrodynamical problem but rather approximate the flow 
evolution by considering that fast material is "accreted" by 
the slow shell. The accretion rate is given by 



^ = M(1 
dt K 



7 ) 



(16) 



where t is the observer time and 7 = T/Ti (F and M being 
the current Lorentz factor and mass of the slow shell). Due 
to the accretion of fast moving material, the Lorentz factor 
of the slow shell increases. When a mass element AM is 
accreted the Lorentz factor becomes 



1/2 



\ TdM + TiM J 

so that 

d7 _ 1 — 7 2 
dM ~ 2M ' 

which can be integrated to give 



1+7 
l- 7 



1 + 7Q 
1 - 7o 



(17) 



(18) 



(19) 



where fj, = M/M and 70 = To/Li. Introducing to = M /M 
and t = t/to, eqs ll6tlT^l vield 



^ = Q(l-7 2 )(l-7) 2 



(20) 



with Q = 



1 ( 1 + 70 N 

2 I 1-70^ 



Equation |201 has the analytical solution 

(21) 

where the function ^(7) is given by 



T= Q [i?(7) " F(7o)] 



1-7 



4(1-7) 4(l-7) 2 
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Figure 4. Pulse profile N(t) (full line) and 1/N(t) (dashed line) 
produced by a relativistic outflow with Ti = 400 decelerated by 
a slow shell of initial Lorcntz factor Tq = 100 (see text for the 
other model parameters). The dashed vertical lines limit the time 
interval where we plot the HIC and the HFC in Fig. 



The solution 7(7") corresponding to ea. l21l has been rep- 
resented in Fig. [5] for 70 = 0.25. When T>2, it is well ap- 
proximated by 

7 (t) * 1 - -±= • (23) 



Once j(t) is known it is possible to calculate the dissi- 
pated power 



f(T) = ^L(l- 7 2 )(l- 7 f 



(24) 



which has been represented in Fig. El At large r, it behaves 
as r -3 / 2 since 

£{t) oc (l- 7 ) 3 (l + 7) oct" 3/2 (25) 

for t>2. Pulse evolution is essentially completed at r ~ 10 
when F/Fi > 0.8 and £ has decreased by more than an 
order of magnitude. 

The dissipated power given by eq. l24l is slightly different 
from what the observer will see since the energy released at 
time t is spread over an interval At corresponding to the 
difference in arrival time for photons emitted by a shell of 
radius r moving at a Lorentz factor V 



At 



2cT 2 



(26) 



The solution for £ including angular spreading has been 
obtained numerically and is also shown in Fig. |^] It differs 
from the analytical expression (eq. 1241 at early times but 
preserves the power law decay of slope e = 3/2 at late times. 



4 SPECTRAL EVOLUTION AND EMISSION 
PROCESSES 

We now use the analytical model to follow the spectral evo- 
lution during pulse decay. If the dissipated energy is radiated 
by the synchrotron process the peak energy E p is 
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Figure 5. Hardness-intensity (full line) and hardness-fluence 
(dashed line) correlations for the pulse shown in Fig. |1] 




(27) 



where B is the magnetic field and F c the characteristic elec- 
tron Lorentz factor behind the shock. With classical equipar- 
tition assumptions B and r c can be expressed as 

2\l/2 



B — (SvraB/oeC 
and 

r. 



a c m p 



(28) 



(29) 



where p is the density and ec 2 the dissipated energy per 
unit mass (both in the comoving frame) ; «b and a a are the 
equipartition parameters and £ is the fraction of electrons 
which are accelerated. Finally, 



1/2 5/2 



E Byn (X T p 

where the comoving density p is proportional to r 
the shock radius r ~ T 2 ct) and e is obtained from £ 
^Tec 2 and eq.[H 



(30) 
" 2 (r being 



(1-7) 2 



2 7 

This leads to the following expression for E syn 

c: q-7) s 



(31) 
(32) 

which behaves as a power law (E p oc t~ 7 ^ 2 ) when (1 — 7) ~ 
£-1/2 This is much steeper than the observed spectral evo- 
lution of pulses which satisfy both the HIC and the HFC. 
Instead of using eq. 1301 we therefore parametrize the peak 
energy with the more general phenomenological expression 

l) 2v 



E p ccTpV oc 
which becomes 



J! 

,y4x+y — l-f-2x 



(33) 



(34) 



at late times. The exponents x and y can be different from 
their standard synchrotron values 1/2 and 5/2 if the equipar- 
tition parameters ob, a e or C vary wi th p or/and e. For ex- 
ample, iDaigne fc Mochkovitchl jl998i) adopted a fraction £ 
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of accelerated electrons proportional to e so that F c remains 
constant which leads to x = y = 1/2 and E p oc t _3//2 . How- 
ever, since most of the observed values of <5* = 2x + y are 
smaller than 1.5, it seems necessary to further reduce the x 
and y indices and we have therefore considered below the 
case x = y = 1/4, i.e. <5* = 0.75. 



5 TEMPORAL PROFILES 

We obtain the temporal profile of synthetic pulses from 
egsl9"l l2"Il and l33l of our model. We have represented in Fig. 2] 
a pulse formed when a wind of Lorentz factor Ti = 400 and 
power MV lC 2 = 10 52 erg.s is decelerated by a slow shell 
with To = 100. We adopt to = 0.4 s, x = y = 1/4 and z = 1. 
The profile is computed in the BATSE range (20 - 1000 keV) 
and the constant of proportionality in eg. 1331 is fixed to get 
a peak energy E p = 300 keV for the whole pulse spectrum. 
The pulse duration is close to 10 (1 + z) to as expected from 
the results obtained in Sect. 3. 

The evolution after maximum is initially close to a 1/t 
decay (i.e. n ~ 1 in eq. |SJ as can be seen in Fig. 2] where 
1/N(t) has been also represented. This can be simply un- 
derstood from eq. 1141 which, for the decay slopes of the 
dissipated power e — 1.5 (eq. 1251 and of the peak energy 
5* = 2x + y — 0.75 gives 



n = -0.75 C 



(35) 



With — 2 < ( < — 1, the central value of n is indeed close 
to unity. Since the decay phase of our synthetic pulse can 
be decribed by eqs. 2] and GD it should also satisfy both the 
HIC and the HFC. This is checked in Fig. [K] where the two 
relations have been plotted in the time interval delimited by 
the two vertical lines in Fig.0] The HIC is not a strict power 
law but its average slope S ~ 0.9 is in global agreement with 
the observations. The HFC is satisfied to a better accuracy 
since the relation between log E p and the photon fiuence is 
quasi-linear in the considered interval. 



6 TIME LAGS 



iNorris et alJ <l2000f) have shown that time lags between dif- 
ferent energy channels correlate with spectral hardness and 
possibly also with the burst peak luminosity. GRBs are dis- 
tributed in a triangular domain of the time lag-hardness 
ratio diagram (the hardest bursts having the shortest time 
lags; see Fig. |SJ and the time lag-luminosity relation ob- 
tained from 6 bursts with known redshifts takes the form 

At 3 i 



130 



01 s 



(36) 



where At3i is the time lag between BATSE channels 3 and 
1 and Z/51 is the luminosity in units of 10 51 erg.s -1 . In our 
model, we estimate time lags by cross-correlating profiles 
in different energy channels obtained from eq. |2l The first 
factor in eq. |5Jis 



E P (t) 



(37) 



which behaves as i 2x +y~ 1 --' during pulse decay. The sign of 
A = 2x + y — 1.5 is of great importance in determining the 
time tmax of maximum count rate and the related time lags. 
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Figur e 6. Time l a g — ha rdness ratio diagram of BATSE bursts 
from INorris et alJ l2000l) compared to model predictions. The 
three lines represent sequences of synthetic pulses obtained with 
similar distributions of the Lorentz factor (the dotted line corre- 
sponds to Vi = 300, the full line to Ti = 400 and the dashed line 
to Ti = 600; see text for details). 



If A < the function tp(t) has a maximum at some early 
time t m before decreasing as t A while it steadily increases 
for A > 0. The second factor in eq. |3]is ip{E p ) so that t ma x 
is solution of the implicit equation 



1>{t 



dlog <p 



dlog E p 



= , 



(38) 



where _E™ ax and i?™ ax are the values of E p and its time 
derivative at t — £ max . Since in most pulses the evolution 
of E p precedes the count rate, i?™ ax is negative. If the low 
energy slope of the Band spectrum a < —1, the derivative 

dlog v is positive (see Fig. ^ and ea. 1381 then shows 

> which, for A < 0, leads to 

(39) 



d log Bp 

that ip(t 

tmax ^ t 



which provides a strict upper limit on the time lags between 
different energy channels 



At < t n 



(If a > -1, * loe £ 



E p but even in this case t 



(40) 

can be weakly negative at large 
nax never greatly exceeds t m .) 



When A > 0, no constr aint such (40) applies and the 
time lags can be quite large. iDaiene fc Mochkovitcbl (|l998h 
who adopted x — y = 1/2 (A = 0) obtained lags of several 
seconds between BATSE channels 1 and 3 (or 4) for a pulse 
lasting about 10 s, while currently observed values are in 
the range 10 -2 to a few 10 -1 s. Moreover, for A > the 
time lags increase with pulse hardness, in contradiction 
with the observations. 

Conversely, with x — y = 1/4 the time lags are short 
(At3i < 0.5 s) even for long pulses and they decrease with 
increasing hardness. Figure |S| shows the time lag - hardness 
ratio relation given by our model superimposed to the 
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At 31 (s) 

Figure 7. Time la g — luminosi t y cor relation predicted by our 
model compared to lNorris et alj l200(j) results for 6 GRBs with 
known redshifts. The lines correspond to the three cases already 
considered in Fig, IE1 

iNorris et all fcOOd) results for BATSE bursts. The thick 
grey lines in Fig. HJ correspond to sequences of pulses of 
comparable duration (tgo ~ 10 s) obtained with similar 
distributions of the Lorentz factor (a slow shell with 
To = 100 decelerating a fast wind with Ti = 300, 400 
or 600) and a varying value of the (isotropic) injected 
power (from 5 10 51 to 10 54 erg.s -1 ). Even if differences in 
duration and redshift will contribute to add scatter to their 
distribution it appears that synthetic pulses populate the 
same triangular domain as observed ones. 

We finally checked if our model was able to reproduce 
the time lag - luminosity correlation (eq. 1361 . When x — 
y — 1/4, we do find that the lags decrease with increasing 
luminosity. This is a consequence of the HIC and the time 
lags - hardness ratio relation discussed above. The results 
are shown in Fig. |7| where the three lines correspond to the 
wind cases with Fi = 300, 400 and 600 already considered 
in Fig. |S| It can be seen that there is an overall agreement 
between the model predictions and eg. 1361 However, at low 
luminosities and large time lags we obtain a rather wide strip 
instead of a single relation as eq. 1361 If this is confirmed by 
the analysis of more GRBs with known redshifts, the time 
lag - luminosity correlation will still be useful for statistical 
studies of large burst samples but may be quite inaccurate 
to estimate the luminosity of a specific event. 



7 CONCLUSION 

We have developed a simple model where GRB pulses 
are produced when a rapid part of a relativistic outflow is 
decelerated by a comparatively slower shell. We do not solve 
the true hydrodynamical problem but rather assume that 
the slow shell "accretes" the fast moving material which 
allows to obtain an analytical solution for the dissipated 
power £ . During pulse decay £ oc t~ 3 ^ 2 as £ oc t~ 3 when 
the evolution is fixed by shell geometry. To compute the 
spectral evolution of our synthetic pulses we parametrize 
the peak energy as E p oc p x e y F where p, ec 2 and T arc 



respectively the post shock values of the density, dissipated 
energy (per unit mass) and Lorentz factor. At late times, 
we get E p oc t~ i - 2x+v ^ which constraints x and y since in 
most observed bursts -Bp oc t~ s * with 5* < 1.5. The syn- 
chrotron process with standard equipartition assumptions 
corresponds to x = 1/2 and y = 5/2 (i.e. 2x + y = 3.5) and 
gives a much too steep spectral evolution. One has then to 
suppose that the equipartion parameters a c , qb and £ vary 
with p or/and e to reduce x and y (a possible alternative 
being that energy is radiated by another process - different 
from synchrotron - but which can still be approximated by 
ea.1331. 

We have considered the case x = y = 1/4 and the re- 
sulting pulses then have temporal and spectral properties 
in excellent agreement with the observations. They follow 
both the HIC and the HFC during the decay phase and the 
time lags between energy channels decrease with increasing 
pulse hardness and peak luminosity. We therefore conclude 
that if GRB pulses are produced by internal shocks, their 
temporal and spectral properties are probably governed by 
the hydrodynamics of the flow rather by the geometry of the 
emitting shells. 



REFERENCES 

Band D., Matteson J., Ford L., Schaefer B., Palmer D., 
Teegarden B., Cline T., Briggs M., Paciesas W., Pendleton 
G., Fishman G., Kouveliotou C, Meegan C, Wilson R., 
Lestrade P., 1993, ApJ, 413, 281 
Daigne F., Mochkovitch R., 1998, MNRAS, 296, 275 
Dermer C. D., Mitman K. E., 1999, ApJ, 513, L5 
Fenimore E. E., Madras C. D., Nayakshin S., 1996, ApJ, 
473, 998 

Fishman G. J., Meegan C. A., Wilson R. B., Brock M. N., 
Horack J. M., Kouveliotou C, Howard S., Paciesas W. S., 
Briggs M. S., Pendleton G. N., Koshut T. M., Mallozzi 
R. S., Stollberg M., Lestrade J. P., 1994, ApJS, 92, 229 
Golenetskii S. V., Mazets E. P., Aptekar R. L., Ilinskii 

V. N., 1983, Nat, 306, 451 
Granot J., Piran T., Sari R., 1999, ApJ, 513, 679 
Kobayashi S., Piran T., Sari R., 1997, ApJ, 490, 92 
Liang E., Kargatis V., 1996, Nat, 381, 49 
Meszaros P., Rees M. J., 1993, ApJ, 405, 278 
Norris J. P., Marani G. F., Bonnell J. T., 2000, ApJ, 534, 
248 

Norris J. P., Nemiroff R. J., Bonnell J. T., Scargle J. D., 
Kouveliotou C, Paciesas W. S., Meegan C. A., Fishman 
G. J., 1996, ApJ, 459, 393 

Rees M. J., Meszaros P., 1994, ApJ, 430, L93 

Ryde F., Svensson R., 2000, ApJ, 529, L13 

Ryde F., Svensson R., 2002, ApJ, 566, 210 

Soderberg A. M., Fenimore E. E., 2001, in Costa E., Fron- 
tera F., Hjorth J., eds, Gamma-ray Bursts in the After- 
glow Era Berlin Heidelberg: Springer, p. 87 



